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A model approach to the description of static stars filled with a charged Pascal perfect fluid within the framework 
of general relativity is investigated. The metric is written in Bondi's radiation coordinates. The gravitational 
equations are reduced to a nonlinear oscillator equation after transfomation to a new variable as a function of the 
radial coordinate. It is shown that in this case exact solutions of the Einstein-Maxwell equations for a concrete 
energy density distribution law of the charged fluid may be obtained as solution of the harmonic oscillator equation. 
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In general relativity, together with Schwarzschild's 
exterior solution (see, e.g., [1]), which describes the 
gravitational static field of spherical objects, there is 
a solution of the gravitational equations for an exte- 
rior field of a massive spherical body with an elec- 
trical charge, well known as the Reissner- Nordstrom 
solution (see, e.g., [2]). However, as to a similar prob- 
lem of finding of the solution for an interior field of an 
electrified body, there are quite clear dimcultes of the 
solution of the Einstein-Maxwell nonlinear differential 
equations in a medium. 

Let us consider an interior static stellar model hav- 
ing some distribution of the electric charge, within the 
framework of general relativity. The exterior gravita- 
tional field of such a star is the Reissner- Nordstrom so- 
lution [2]. The metric is written in Bondi's radiation 
coordinates: 



ds 2 = F(r)dt 2 + 2L{r)dtdr ~ r 2 (d6 2 + sin6 2 d(p 2 ), (1) 

where t is a time coordinate, r, is a radial coordinate, 8, 
and ip are angular variables; the Greek indices take the 
values 0,1,2,3; F(r) and L(r) are metric functions 
to be found. They depend only on the radial variable. 
The velocity of light and the Newtonian gravitational 
constant are chosen to be equal to unity. 

The Einstein equations for the interior fielf are writ- 
ten with a source for a charged substance as 



G a p = R a p — -j9apR = — aeT a ^, 



(2) 



where G a p is the Einstein tensor, R a p is the Ricci ten- 
sor, R — R a a is the scalar curvature, k = 87r is Ein- 
stein's constant, g a p is the metric tensor. The resulting 
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energy-momentum tensor of the charged substance T a p 
is taken here as a sum of the energymomentum tensor 
of the Pascal neutral perfect fluid 
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fluid 
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and the electromagnetic energy-momentum tensor 
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= —(-F a ,F»+g aP F^Fn, 
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where /i(r) is the density, p is the pressure, u a — 
dx a /ds is the 4- velocity, F a p is the electromagnetic 
field tensor, all functions depend only on the radial vari- 
able. 

The Maxwell equations (the second pair) in our case 
become 
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where g is the determinant of the metric tensor, is 
the electrical current density. 

The statics nature of the system allows one to intro- 
duce a fixed comoving frame of reference with 
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and the physically observable quantities: the mass- 
energy density, the electrical charge density, the elec- 
tric intensity (radial component) and the electric field 
energy density will be written accordingly as 



Uphys = MM; P P hys = P = fu^= j'°v / 5oo; (7.1) 
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As a result, the set of Einstein-Maxwell equations 
is reduced to the following set of nonlinear differential 
equations: 



F 



(InL)' = ae(p + /z); 



F" + —F' - F'{lnL)' 



= -x(p + W el ); 



(8.1) 
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+ W e / ; (8.3) 
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where the prime denotes a derivative with respect to 
x = r/R, < x < 1, R is the radius of the star, 
VC e i = E 2 /8ttL 2 , E is the electric field intensity. 

By the replacement e = F/L 2 it is possible to derive 
a linear second-order ordinary differential equation with 
variable coefficients from the first three Einstein equa- 
tions by eliminating the presure and the mass density: 



G" + f(x)C'+g(x)G = 0, 



(9) 



where f(x) = (Imp)' , ip = y/e/x, g(x) = [2(1 — e) + 
xe']/(2x 2 e)-2 x W e i/e; X = *R 2 - 

For the new variable Q = ((x) introduced by the 
relation 

dC = -^= = -^U (io) 

Eq.(9) is rewritten in a self-conjugate form and becomes 
the equation for a nonlinear oscillator 



G'l ( + n 2 (ax))G = o, 



(11) 



where y 



ft 2 is the squared " frequency" . 



Eq.(10) not always can be integrated in elementary 
functions, and more simply the function £l 2 is written 
as a function of the variable y : 



= -d(<S>/y)dy-2 X W el /y. 



(12) 



Taking into account the contribution of the electric 
field energy, the function $ playing the role of the New- 
tonian gravitational potential inside the star can be ob- 
tained from the gravitational equations as 

$ = 1 - E = | J (H( X ) + W e l )x 2 dx 

= ^=f(»(y) + w el )Vydy. (13) 



Choosing the mass density /i and the electric field 
energy density inside the star W e i as square-law func- 
tions of x (as was made in [3]), 



H = /u (l - x 2 ) = no(l - y): 
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(14) 
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we obtain from (11) the equation for a linear oscillator 

G' c ' c + n 2 G = (16) 
with a constant value of the natural frequency 

(17) 



o2 X, HA 2 



5 vr " 8tt 

The continuity of the electric energy density on the 
stellar surface requires A = Q/R 2 , where Q is the 
integral electric charge of the star and /j,q is its central 
mass density. 

In this case the function $ may be written as 



,/i 1/A 2 , , 



(18) 



The radial dependence of the electric charge density 
in the interior field of the star can be found for the 
chosen functional law of the observable electrical energy 
density from Maxwell's equation (8.4): 



P (x) = Q ■ yfi(x)/V, (19) 

where V = 47ri? 3 /3 is the Euclidean volume of 3-space. 

The general solution of Eq.(16) for the required func- 
tion G will be written as 

G = G cos(O C + a), (20) 
and for the metric function g ao = F = G 2 

F = G 2 cos 2 (n ((x)+a), (21) 

where ((x) is a quadrature depending on the relation 
between the parameters and A. 
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